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Abstract 

Large high-quality single crystals of Bi2Sr2Ca n _iCu n 02n+4+<5 in 
= 2, 3) were successfully grown using an improved traveling solvent 
floating zone (TSFZ) method, which features a slow growth rate and 
steep temperature gradient along the melting zone. By measuring 
anisotropic resistivities and susceptibilities of Bi2Sr2CaCu20s+<5, the 
characteristic pseudogap temperature T* was studied as a function 
of doping. The T* suggest that the pseudogap is not simply a pre- 
cursor of high-T c superconductivity, but that the pseudogap and the 
superconducting gap compete with each other. The anisotropic resis- 
tivities of Bi2Sr2Ca2Cu30io+<5 were also measured, revealing that the 
T c remains fixed in the over doped region while anisotropy decreases 
continuously. This anomalous behavior will be discussed in terms of 
the inequivalent hole doping, which occurs between two inequivalent 
CuC>2 planes in the triple-layer system. 
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1 INTRODUCTION 



Ever since high-T c superconductivity was discovered more than 15 years ago 
PQ, intensive investigations have been carried out in attempts to reveal its 
mechanism. However, the mechanism of high-T c superconductivity has re- 
mained elusive because we have not been able to resolve some important 
issues in its physics. Resolving these issues has been difficult because of the 
lack of large high-quality single crystals, which are needed so that we can 
examine the anisotropic properties more precisely. In a case that the crys- 
tal is a new one, like Bi 2 Sr2Ca 2 Cu30io+«5 (Bi-2223) here, the benefits of the 
successful crystal growth is especially large. 

One unresolved issue in the physics of high-T c superconductivity is the 
pseudogap that depletes the electronic density of states (DOS) around the 
Fermi level from characteristic temperatures (T*) far above the supercon- 
ducting transition temperature (T c ). Since the pseudogap evolution is com- 
monly observed in most high-T c superconductors [2|, it is generally believed 
that the pseudogap is closely related to the pairing mechanism of high-T c su- 
perconductivity. Therefore, the pseudogap phenomenon has been attracting 
much attension and a huge amount of experimental data has been accumu- 
lated. However, up to now, as shown in Fig. 1, there is still no consensus 
on when (at what temperatures) the pseudogap opens. One of the reasons 
for the lack of consensus is that there is probe dependence. Many probes 
such as angle resolved photoemission spectroscopy (ARPES) Ej, NMR 
IE] and in-plane resistivity p a [Zj suggest that the T* is close to line "A". 
On the contrary, the uniform susceptibility % °f La 2 _ a; Sr 3 ;Cu04 [Hj suggests 
that the T* is close to line "B". Another reason is that there is material 
dependence. For example, the x °f YBa 2 Cu 3 6+( 5 ;9j suggests that the T* is 
close to line "A". Clarifying whether the pseudogap phase boundary T* is 
"A" or "B" is crucially important because it would elucidate the relationship 
between the pseudogap and superconducting gap. The T* of boundary "A" 
implies that the pseudogap is a precursor of high-T c superconductivity, while 
that of the boundary "B" implies that the pseudogap and superconducting 
gap compete with each other. To find the generic pseudogap phase bound- 
ary, the Bi 2 Sr 2 CaCu 2 0s+,5 system has been studied in detail using its large 
high-quality single crystals. 

Another unresolved issue is the physical origin of the correlation between 
the maximum T c and the number (n) of Cu0 2 planes. Empirically, the max- 
imum T c first increases with n for n < 3, then decreases JU]- Taking the Bi 
family, i.e. Bi 2 Sr 2 Ca n _iCu n 2n+ 4 + <5 (n = 1-3), as an example, the maximum 
T c is approximately 34, 90, and 110 K for n = 1 (Bi 2 Sr 2 CuO6+<5:Bi-2201), n 
= 2 (Bi 2 Sr 2 CaCu 2 8+5 :Bi-2212), and n = 3 (Bi 2 Sr 2 Ca 2 Cu 3 Oi 0+ 5:Bi-2223), 
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respectively. Single crystals of triple-layered Bi-2223, however, had been 
very difficult to grow and large samples for physical property measurements 
were therefore unavailable for a long time. The difficulty is attributed to the 
narrowness, both in temperature and composition, of the liquidus line for 
Bi-2223. To overcome this serious problem, we have developed an improved 
traveling solvent floating zone (I-TSFZ) method. 

The next section describes the concept of the I-TSFZ method and its suc- 
cessful application to Bi 2 Sr 2 Ca n _iCu ri 02 n +4+<5 (n = 2, 3) systems. Section 3 
addresses, based on the systematic measurements of those crystals, several 
important issues on high-T c superconductivity: the pseudogap phase bound- 
ary, the origin of semiconducting out-of-plane transport p c (T), the origin of 
anomalous temperature/doping dependence of magnetic uniform susceptibil- 
ity Xi an d the authenticity of charge confinement. In addition, for Bi-2223, 
anomalous doping dependence, in which T c does not decrease with increas- 
ing doping in the overdoped region, will be revealed. The reason for the 
anomalous behavior as well as the reason why T c is higher in triple-layered 
compounds will be discussed. Section 4 summarizes the paper. 

2 CRYSTAL GROWTH 

2.1 FLUX METHOD vs. TSFZ METHOD 

Growing large high-quality single crystals is very important for reliable phys- 
ical property measurements. Since all high-T c superconductors melt incon- 
gruently, i.e. decompose before melting, direct crystallization from the melt 
(melt growth), such as occurs in the Bridgman method or the Czochralski 
method, can not be used because it results in some other phase. The al- 
ternative is solution growth, in which the constituents of the material to be 
crystallized are dissolved in a suitable solvent and crystallization occurs as 
the solution becomes critically supersaturated ^Tj. Among solution growth 
methods, the most popular two in the crystal growth of the high-T c supercon- 
ductors are the flux method (slow cooling method) and the traveling solvent 
floating zone (TSFZ) method. 

The principle of the flux method is shown in Fig. 2 as (a) a phase diagram 
and (b) a practical image. In this method, the supersaturation is obtained by 
cooling. Continuous crystal growth occurs because the solubility decreases 
with decreasing temperature as shown in Fig. 2(a). The flux method is 
simple and inexpensive and is thus widely used in preparing new phases. On 
the other hand, the grown crystals are often contaminated from the crucible 
used and it is difficult to grow large crystals. More seriously, crystals cannot 
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practically be obtained in a system where the solubility difference available 
is very small (i.e., the liquidus line is very steep and/or narrow). 

The principle of the TSFZ method is shown in Fig. 3 phase 
diagram and (b) a practical image. A proper amount of solvent with the 
composition "s" is set between a feed rod and a seed crystal, melted, and 
then they are slowly moved downward [Fig. 3(b)]. Consequently, crystals 
with the composition of "T" are precipitated on the seed crystal through 
a solid-liquid reaction between "T" and "s". The solvent composition "s" 
is maintained unchanged because the solute exhausted by the crystalliza- 
tion is continuously supplied (transported) from the feed by melting in the 
steady growth state. The supersaturation, in this case, is obtained by nutri- 
ent (solute) transport. Therefore, crystal growth continues at one point on 
the temperature-composition phase diagram as shown in Fig. 3(a). Conse- 
quently, the TSFZ method allows us to grow large single crystals even if the 
liquidus line (crystallization field) is very narrow. For the same reason, this 
method is especially useful for the preparation of homogeneous solid solutions 
and homogeneously doped crystals [12]. Furthermore, the crystal growth is 
performed in a free space, i.e., a crucible is not used, so the obtained crystals 
are free from impurities. It is difficult, however, to keep the molten zone sta- 
ble, because it is supported only by surface tension. For example, the molten 
zone for various reasons, sometimes drops down, soaks into the feed rod, and 
bubbles by many reasons. Therefore, the applicability of the TSFZ method 
has generally been thought to be narrow. We choose the TSFZ method for 
the crystal growth of Bi 2 Sr2Ca n _iCu n 02n+4+5 (n = 2, 3) systems, because, 
in principle, it makes growth of the large high-quality single crystal possible. 

2.2 CONCEPT OF THE I-TSFZ METHOD 

As mentioned above, some materials are difficult to obtain in single crystal 
form, mainly because they have only narrow liquidus lines. With narrow 
liquidus lines, the supersaturation is small, which makes it difficult to grow 
large crystals by standard crystal growth methods such as the flux method. 
Figure 4 shows a schematic phase diagram of the Bi family. Obviously, 
Bi 2 Sr 2 Ca n _iCu n 02n+4+<5 (n = 2, 3) systems, especially the n = 3 compound, 
Bi-2223, has a narrow liquidus line. To overcome this serious problem of 
the insufficient supersaturation, we developed the I-TSFZ method with two 
features [HUE]. 

One feature of the I-TSFZ method is the very slow growth rate adopted. 
The driving force for crystal growth is provided by supersaturation ^T] . The 
maximum growth rate increases as the level of supersaturation increases. 
Therefore, if the supersaturation level is very low, we should use a very 
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slow growth rate so as not to exceed the maximum growth rate. Further- 
more, these Bi compounds (Bi-2212 and Bi-2223) are known to grow very 
anisotropically; i.e., the growth rate along the c-axis is extremely low com- 
pared with that along the a- or 6-axis. Probably because of this anisotropic 
growth rate, single crystals of these compounds are usually thin (< 0.1 mm 
and < 0.003 mm for Bi-2212 and Bi-2223, respectively) along the c-axis. The 
very slow growth rate is expected to allow sufficient time for growth along 
the c-axis and thus increase the crystal thickness. 

The other feature is the steep temperature gradient at the solid-liquid 
interface. In the TSFZ method, the whole temperature range covering the 
liquidus line is not necessarily available because the molten zone sometimes 
drops depending on its size, viscosity, and other factors. We found that a 
steep temperature gradient provides comparatively flattened solid-liquid in- 
terfaces for both feed and grown crystal sides, which makes it possible to 
reduce the length of the molten zone. It is considered that the molten zone 
becomes stable if its size is reduced because it is supported only by surface 
tension. Therefore, the steep temperature gradient is expected to stabilize 
the molten zone and thereby make it possible to control the growth temper- 
ature precisely. On the other hand, it is important in the TSFZ method that 
the solvent composition and quantity be kept exactly at the specific optimal 
values. Fortunately, in the TSFZ growth of Bi-based superconductors, the 
solvent composition can be self-adjusted during growth if the molten zone 
is sufficiently stable against composition and temperature fluctuation. The 
steep temperature gradient is expected to enhance the self-adjustability by 
stabilizing the molten zone. In addition, more importantly, the supersatura- 
tion level increases with increasing temperature gradient. As noted above, 
supersaturation is the driving force of crystal growth in solution growth. 
Therefore, the steep temperature gradient is expected to enhance the maxi- 
mum growth rate, which was considered to be relatively lower in the system 
having narrow liquidus line (crystallization field). 

Noteworthy is the intimate relation between the conditions for stable 
crystal growth and the I-TSFZ method. The instability of a planar crystal 
surface growing in a doped melt or in a solution has often been discussed 
in terms of constitutional supercooling ^3]. As the crystal grows, solvent is 
rejected at the crystal surface and so the solute concentration immediately 
ahead of the interface becomes appreciably lower than that in the bulk of the 
solution [Fig. 5(a)]. This accumulation of solvent near the interface results in 
a depression of the equiliblium liquidus temperature T L as shown in Fig. 5(b). 
If the actual temperature distribution in the solution is as shown in line "A" 
of Fig. 5(b), any protuberance on the interface will tend to grow since it will 
experience higher supercooling. The region ahead of the interface will then 
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become unstable, which results in the so-called cellular growth with multi- 
domained crystals and precipitations in between them. Stable growth can 
occur if the temperature gradient at the interface is increased as in " B" of Fig. 
5(b), in which case the actual temperature everywhere ahead of the interface 
is higher than the liquidus temperature. The condition for constitutional 
supercooling when it is applied to the solution growth is expressed by Tiller 
PS| as 



where v is the growth rate, D a the diffusion coefficient of the solvent, and 
rrii, k*, n iy and Dj are the slope of the liquidus curve, the effective partition 
coefficient, the concentration in the solution, and the diffusion coefficient of 
the solute constituent i, respectively. The effective partition coefficient is 
defined by 



where n l c and n\ n are the concentrations of each solute constituent i in the 
crystal and solution, respectively. From eq. (1), the very slow growth rate 
and the steep temperature gradient adopted in the I-TSFZ method is ex- 
pected to work cooperatively to prevent the occurrence of constitutional su- 
percooling, which is a cause of cellular growth. 



Single crystal growth of the Bi-2212 phase by the TSFZ method was first 
reported by Takekawa et al. . For the crystal growth, they used an 
infrared convergence-type floating zone furnace with a 1500-W halogen lamp 
as the radiation source. The growth rate was 1 mm/h. The grown boule 
was an aggregate of large single crystal grains with growth direction along 
the Although the size was sufficient in a- and 6-axis direction (several 

mm), that in the c-axis direction (thickness) was limited below 0.1 mm. 

To increase the crystal thickness, we decreased the growth rate to 0.5 
mm/h. In addition, we increased the temperature gradient by using 300-W 
halogen lamps (A furnace equipped with two halogen lamps was used) in 
place of the standard 1500-W lamps. The temperature gradient was about 
370°C/cm, which is two times larger than that with 1500-W lamps. Actu- 
ally, the length of the molten zone is reduced from 4.0-4.5 mm for 1500-W 
halogen lamps to 3.5-4.0 mm. Compositions of Bi 2 .iSr! . 9 CaCu 2 8 + ( 5 and 
Bi2.2Sri. 6 Cao.85Cu 2 .20 8+( 5 were used for the feed and the molten zone, respec- 
tively. The samples were prepared from powders of Bi 2 03, SrCOs, CaCOs, 




(1) 




(2) 



2.3 Bi 2 Sr 2 CaCu 2 8+( 5:Bi-2212 
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and CuO (all of 99.999% purity) mixed in the desired ratio, and then calcined 
at 780°C (feed) and 760°C (molten zone) for 12 h with repeated regrindings. 
The calcined powders were again reground, placed into a rubber tube and 
then hydrostatically pressed under 2000 kg/cm 2 . The pressed rods for the 
feed and for the molten zone were sintered at 860°C and 830°C, respectively, 
for 50 h. The sintered feed rod was premelted at a rate of 25 mm/h to pre- 
pare a dense feed rod. The feed rod was suspended at the bottom of the 
upper shaft, and a seed crystal (or just a sintered rod) was held at the top of 
the lower shaft. A proper amount of sintered rod (pellet) for the molten zone 
was set under the bottom of the feed rod and was melted to form a proper 
solid-liquid reaction. The upper and the lower shafts were counterrotated at 
a rate of 10.5 and 10.0 rpm, respectively. Then, the lower shaft was moved 
downward at a rate (growth rate) of 0.5 mm/h, while the upper shaft was 
moved at a rate of 0.4 mm/h. The growth atmosphere was a mixed gas flow 
of 2 (20%) and Ar (80%). 

Figure 6(a) shows an image of the whole grown boule. Crystals tend to 
become large as growth proceeds. Figure 6(b) shows an optically polarized 
micrograph of a cross section normal to the growth direction of nearly the 
last part of the grown boule. As can be seen in this figure, large crystals with 
the thickness of about 1 mm were obtained. Figure 6(c) shows an afr-plane 
image of such crystals cleaved from the boule. The afr-plane of the crystal 
was fairly large: the length along a-axis was over 5 mm and that along fe-axis 
was 4.5 mm, which is almost the thickness of the grown boule. The main 
reasons for this successful result are that enough time was allowed for growth 
along the c-axis by the slow growth rate and that the c-axis maximum growth 
rate was enhanced by steep temperature gradient. These crystals have been 
used for the physical property measurements (section 3). A further increase 
of the temperature gradient, unfortunately, caused unstable crystal growth. 
Prior to our result, a similar or better result of 1.5 mm thick Bi-2212 crystal 
had been reported by Gu et al. [IB], who used a growth rate of 0.2 mm/h 
and a halogen lamp power of 500 W. 

In an attempt to further enlarge the crystals, the growth rate was re- 
duced to 0.05 mm/h with the other conditions fixed \F3\. Figure 7(a) shows 
a cross-sectional image of the grown boule obtained with an optically po- 
larized microscope. About one-third of the crystal boule (left part) was a 
single domain with thickness of 1.5 mm, while the remaining part 

was heavily multi-domain. The obtained Bi-2212 crystals unfortunately con- 
tained a large number of Bi-2223 impurities. That such thick crystals were 
grown indicates that the slow growth rate is indeed effective for enlarging 
crystal size. However, the existence of heavily multi-domain region, on the 
other hand, suggests that another factor competes with crystal enlargement. 
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Reflecting the extreme anisotropy in the growth rate between the b- axis and 
c-axis, the cross section becomes elliptical in shape as shown in the figure. 
This ellipsoidal cross-section of a crystal boule likely causes the molten zone 
to drop down along the minor axis of the ellipsoid. The situation is illustrated 
in Fig. 7(b). To prevent the molten zone dropping down, the remaining part 
[right part in Fig. 7(a)] may have self-assembled into a heavily multi-domain 
structure. Indeed, the solid-liquid interface of the remaining part is rather 
flat compared to that of the single-domain part [Fig. 7(c)]. Figure 7(d) and 
(e) show enlarged images of the solid-liquid interface of the single-domain 
part and that of the multi-domain part, respectively. A smooth interface is 
seen in Fig. 7(d), while a rough one is seen in Fig. 7(e), evidencing the single- 
domain growth and multi-domain growth, respectively. It should be noted 
that this multi-domain growth is not the cellular growth caused by consti- 
tutional supercooling, which was discussed in section 2.2. Instead, it may be 
due to molten zone instability caused by the anisotropic growth rate. In this 
context, an interesting attempt has been reported by Watauchi et al. [20], i n 
which the anisotropic growth rate was controlled using an infrared heating 
furnace equipped with four mirrors. Taking CuGe0 3 and Sr 14 Cu240 4 i as ex- 
amples, they showed that the ellipsoidal cross-section of a cylindrical crystal 
grown usually under an isotropic lamp power changed to a circular one when 
the crystals were grown under an anisotropic lamp power. 

2.4 Bi 2 Sr2Ca2Cu 3 Oio + 5:Bi-2223 

In initial attempts at growing Bi-2223 phase, the KC1 flux method was used 
[2HI22], but the obtained crystals were very thin, with sizes of 1 x 1 x 0.003 
mm 3 . Obviously, this is due to the narrowness of the liquidus line. To improve 
the size and quality of this system, we adopted the I-TSFZ method [T3|. As 
noted in the section 2.1, the TSFZ method is suitable for such single-crystal 
growth because the growth can be performed continuously at one point on 
the temperature-composition phase diagram. 

Crystal growth of the Bi-2223 phase by the I-TSFZ method was carried 
out with a growth rate of 0.05 mm/h and temperature gradient of 370°C/cm 
using 300-W halogen lamps. For a feed rod, the nominal composition of 
Bi2.iSri. 9 Ca2Cu30io+,5 was used. A solvent was not used. Other growth 
conditions were the same as for Bi-2212 (section 2.3). Figure 8(a) shows a 
feed rod premelted at a rate of 25 mm/h. Making this kind of dense feed 
rod with uniform thickness is indispensable for the successful crystal growth. 
Photographs [Fig. 8(b)] show actual crystal growth near the molten zone. 
They were taken once every week. The crystal growth was very stable for a 
long period (over two months). The final rod is shown in Fig. 8(c). From 
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the last 1-cm part of the grown rod, large high-quality single crystals were 
harvested. Figure 8(d) shows an optically polarized micrograph of a cross 
section normal to the growth direction of this part. Although the grown 
boule consisted of many largely misoriented plate-like crystals of Bi-2223 as 
in the case of Bi-2212, large plate-like single crystals with dimensions of up 
to 4 x 2 x 0.1 mm 3 [Fig. 8(e)] could be cleaved from the upper right side [Fig. 
8(d)]. Figure 9 shows the temperature dependence of the DC susceptibility 
of a crystal also cleaved from the last 1-cm part of the boule. There is a 
sharp superconducting transition (< 3 K) with T c = 105 K. After the sample 
was annealed in flowing oxygen at 600°C, T c increased to the maximum value 
of 110 K, suggesting that the as-grown sample is in the underdoped region. 
The small step-like change in the shielding signal, AM, around T = 83 K 
may have been caused by the intergrowth of the Bi-2212 phase. However, 
AM is less than 2% of the full volume of the demagnetization measured 
at 5 K. This assures that the obtained crystal is a nearly pure-phase Bi- 
2223 with less than 2% Bi-2212 impurity. A detailed structure analysis of 
this sample based on a four-circle x-ray diffraction measurement has recently 
been reported [2Sj- Single crystals in other parts (other than the last 1-cm 
part) were small (typically 1 x 0.5 x 0.02 mm 3 ). They consisted mainly of 
the Bi-2223 phase with around a 20% Bi-2212 impurity phase. Since we 
did not use a solvent, the composition for the molten zone might not have 
been sufficiently suitable at the first to middle stage of the crystal growth. 
However, in the TSFZ growth of Bi-based superconductors, the composition 
for the molten zone is self-adjusted during growth. In the last 1-cm part 
of the grown boule, the composition might have been well stabilized at the 
optimal value for Bi-2223 growth. 

We have successfully grown, for the first time, large high-quality single 
crystals of Bi-2223 by using the I-TSFZ method. A very slow growth rate 
and a steep temperature gradient at the solid-liquid interface made it possible 
to stabilize the growth conditions in the narrow crystallization field of Bi- 
2223. Recently, our result has been reproduced by several groups |2UES]- 
Liang et al. [21] have emphasized that the growth rate affects the phase 
formation of Bi-2223 crystals. They have shown that the as-grown crystals 
obtained at a slow growth rate of 0.04 mm/h consist of a >90% Bi-2223 
phase, whereas as-grown crystals obtained at rates of 0.10 and 0.20 mm/h 
contain predominantly Bi-2212, Ca2Cu03, and a small quantity of the Bi- 
2223 phase. 
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3 SOME PHYSICAL PROPERTIES 



3.1 Bi 2 Sr 2 CaCu 2 8+ j:Bi-2212 

3.1.1 PRECISE OXYGEN CONTROL 

Bi-2212 is a suitable system for studying the electronic phase diagram noted 
in the introduction, because its doping level can be controlled over a wide 
range by varying the oxygen content. However, it had been somewhat difficult 
to control the doping levels, especially into the underdoped region, and retain 
good sample quality To overcome this difficulty, we developed a technique 
to control the oxygen content (5) of Bi-2212 crystals precisely |26j . 

Figure 10 shows an equilibrium phase diagram of Bi-2212 single crystals 
for various 5, which was obtained by a thermogravimetric measurement [25] . 
According to the phase diagram, the oxygen content (5) of the crystals was 
controlled to the required values: the crystals were annealed in a quartz tube 
furnace under the corresponding oxygen partial pressures at 600°C for 10 h, 
slowly cooled («0.4 °C /min) while maintaining the equilibrium oxygen pres- 
sures as a function of temperature to 300-400°C, and then rapidly cooled. 
This process ensured both uniform oxygen content and minimum disorder 
in the crystal. Figure 11 shows several DC susceptibilities of the annealed 
samples as well as that of the as-grown one. The annealed samples show nar- 
rower superconducting transition (<3 K) than the as-grown one (transition 
width «5 K), confirming the uniform oxygen content. The superconducting 
transition temperature (T c ) increased, reached a maximum, then decreased 
with increasing oxygen content (5). The T c 's are plotted as a function of 5 
in Fig. 12. Fig. 13 shows the experimental setup for the annealing. For 
the control of the oxygen partial pressures, a mixed gas flow of O2 and Ar 
at a total amount of 5 1/min was used in the higher oxygen pressure region 
(P02 > 10~ 4 atm). In the lower oxygen pressure region (P G 2 <10~ 4 atm), 
the pressure around the sample was controlled by tuning the conductance of 
the evacuation path using a slotted valve, as well as by tuning the oxygen 
flow while continuously evacuating the quartz tube. A personal computer 
controlled the whole system, i.e. both the oxygen partial pressures and the 
temperatures. Heavily overdoped samples with T c « 60 K (5 ~ 0.3) were 
made separately by high O2 pressure (400 atm) annealing at 500 °C for 50 h 
using a hot isostatic pressing (HIP) furnace. 

3.1.2 p a , p c and x 

Figure 14 shows the temperature dependence of in-plane resistivity p a for 
Bi-2212 single crystals with various oxygen contents (5) [2H]- The overall 
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slope dp a /dT decreases with increasing 5. In a simplest model, resistivity 
p(T) is expressed as 

p( t ) = ( 3 ) 
ne z T{l ) 

where m* is the effective mass, n the carrier concentration, e the charge, 
and t the carrier lifetime. The monotonic decrease of dp a /dT assures that 
the carrier concentrations are indeed increased by the oxygenation process. 
A typical T-linear behavior and a slightly upward curvature of p a are also 
seen in an optimally doped (5=0.255, T c =89 K) and an overdoped sample 
(5=0.28, T c =81 K), respectively. For the underdoped samples (5<0.25), p a 
deviates from high-temperature T-linear behavior at a characteristic tem- 
perature T* a (shown by the arrow in Fig. 14), and decreases rapidly with 
decreasing temperature. In the YBa 2 Cu 3 6+( 5 system, this behavior has been 
interpreted as being due to a decrease in spin scattering, which is caused by 
the pseudogap opening j2BJ- T* a in Bi-2212 is estimated as 220, 200, and 175 
K for 5 = 0.2135, 0.217, and 0.22, respectively. 

Figure 15(a) shows the temperature dependence of out-of-plane resistiv- 
ity p c for a Bi-2212 single crystal (sample I) with various oxygen contents 
[26J. A typical semiconductive temperature dependence is observed in a wide 
doping range. It should be noted that this semiconducting behavior of p c is 
quite anomalous in that the temperature dependence of p a and p c is totally 
different. In general, in a metal, the temperature dependence of p a and p c 
is fundamentally the same (p a (T)/p c (T) = constant) even if the system is 
two-dimensional (2D), whereas their magnitude differs (p a (T)/p c (T) ^ 1), 
reflecting the anisotropic effective mass m*. Although many models have 
been proposed to explain semiconductive p c , there was no consensus [29] . 
Figure 15(b) shows p c for another Bi-2212 single crystal (sample II), focus- 
ing on its optimum and overdoped behavior |2Z|- The p c in an optimally 
doped sample (5=0.25, T c =89 K) shows semiconductive behavior in all tem- 
perature regions measured. With increasing 5, the overall magnitude of p c 
decreases and the characteristic temperature T* c [shown by the arrow in Fig. 
15(b)] for the onset of upturn appears. T* c is estimated as 185, 165, and 130 
K for 5=0.26, 0.27, and 0.28, respectively. For the heavily overdoped sample 
(5~0.3, T c =60 K), p c shows no upturn behavior, suggesting the pseudogap 
completely vanishes. 

Figures 16(a) and (b) show the temperature dependence of magnetic sus- 
ceptibilities Xc (H\\c) and Xab (HJ-c), respectively, of the B-2212 single crys- 
tals with various oxygen contents |2Z|. The overall magnitude of Xc and Xab 
increases with increasing 5. It is presumed that the magnetic susceptibilities 
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can be expressed as 

n 

Xa(T) = Xcore + xY + 9 (ifaXspinih T) (4) 

8=1 

where a denotes c or ab, Xcore denotes the sum of the core diamagnetism 
for the ions present, Xa V * s the anisotropic Van Vleck paramagnetism (or- 
bital contribution from the Cu 2+ ion), Xs P in(i,T) the i-th component of the 
spin susceptibility for the doped Cu0 2 plane (For generality, we assumed n- 
components for the spin susceptibility), and g(i) a denotes their anisotropic 
^-factors. The temperature and doping dependence of \ mainly come from 
the spin susceptibility The susceptibilities in the underdoped (5=0.22, T c =82 
K) and optimally doped (5=0.25, T c =89 K) samples monotonically decreases 
with decreasing temperature. For the slightly overdoped sample (5=0.26, 
T c =86 K), they are almost constant at high temperatures and gradually 
decrease from a characteristic temperature T* [shown by the arrow in Fig. 
16(a)] upon cooling. A further increase in 5 causes a shift of T* to lower tem- 
peratures, whereas, above T* the susceptibilities prominently decrease with 
increasing temperature. T* is estimated as 200, 175, and 135 K for 5=0.26, 
0.27, and 0.28, respectively. The difference in the absolute value between Xc 
and Xab m ay be caused by the anisotropic Van Vleck paramagnetism, which 
is almost completely T-independent jHD]- 

3.1.3 ESTIMATION OF DOPING LEVEL FROM p c (T) 

The above systematic measurement for p c {T) presented us with a new method 
for estimating the doping levels p of Bi-2212. Bi-2212 is now a key mate- 
rial for studying the physics of high-T c superconductivity, since several im- 
portant spectroscopic measurements, such as angle resolved photoemission 
spectroscopy (ARPES) [HUH!, scanning tunneling spectroscopy (STS) [HTJIHSj 
and interlayer tunneling spectroscopy (ITS) [HHUM], can be performed on this 
material owing to its clean surface and highly electronic 2D nature. Deter- 
mining the doping levels (p) of a sample is, therefore, an important issue. In 
general, p is estimated from T c using the empirical relation jSS] 

7f^— = l-82.6(p-0.16) 2 (5) 

-L c,max 

where T CjTnax is the maximum T c value of the Bi-2212 sample. However, in 
this case, controlling the p in a wide range including T C)max is a prerequisite, 
but is not always possible. Alternatively, we found that p c {T) provides a 
simple and reliable estimation for p. 
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Figure 17 plots, for several samples (I, II, and III), ratios of peak value, 
Pc, P eak, of Pc just above T c to p c at room temperature, p c (RT), together 
with their T c , as a function of doping levels (p). Here, excess oxygen S was 
converted to p using the above empirical relation. Although the absolute 
value of p c (T) scatters from sample to sample, the ratio r = p c , P eak / Pc(RT) 
scales in one curve 

0.993 

p = 0.024 + . 6 

F r + 4.696 v ; 

This scaling function enables us to estimate p from a simple p c {T) measure- 
ment. The advantages of this method are that it is possible to estimate p 
using one sample and the estimation is independent of the sample quality. 
This method was actually used in an ITS measurement 



3.1.4 PSEUDOGAP PHASE BOUNDARY 

As noted in the introduction, there is still no consensus on when (at what 
temperatures) the pseudogap opens. This is mainly because of our incom- 
plete understanding of the T dependence of the uniform susceptibility x- 
For example, in the La2- a Sr I .Cu04 system, it is known that x first increases, 
takes a broad maximun, and then decreases from some characteristic tem- 
perature (T*) with decreasing temperature and T* decreases with increasing 
doping (Sr content) [Hj. Our observation for Bi-2212 (Fig. 16) is qualitatively 
similar. D. C. Johnston interpreted the x °f the La2- x Sr x Cu04 system as 
the sum of two-dimensional (2D) 5=1/2 square-lattice Heisenberg antifer- 
romagnetism and T- independent Pauli paramagnetism [S]. Later, there was 
proposed a model based on the mode-mode coupling theory of spin fluctua- 
tions in 2D metals with a technically nested Fermi surface, which explains the 
T dependence as a crossover from the localized to itinerant spin fluctuations 
In those pictures, the gradual decrease in x is not attributed to the for- 
mation of the pseudogap. On the other hand, J. W. Loram et al. have shown 
from susceptibility and high-resolution specific heat measurements that the 
pseudogap develops below T* in the YBa 2 Cu 3 6+ 5 [9\ and La 2 _a;Sr x Cu04 
systems |HZj. Similar results have been reported from NMR Knight shift 
measurements jHEj- Up to now, we have not had a unified description of the 
T/doping dependence of \i that is, that the behavior is pseudogap-like at 
low dopings and Curie-like at high dopings. 

Figure 18 plots the above obtained T* a , T* c and T* together with sev- 
eral other characteristic temperatures, including T c , as a function of carrier 
concentration p |27j . Here again p was estimated from T c using the empir- 
ical equation. T* and T* c coincide for each p. This strongly suggests that 
they have the same origin. Unlike YBa 2 Cu30g+5, which may change from 
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2D to 3D upon doping, Bi-2212 stays in 2D for all doping levels. We have 
previously shown that the p c in the Bi-2212 system is governed by a simple 
tunneling process [12]. Then, p c should reflect the in-plane density of states 
(DOS) through the equation |HH] 



" C tir oo N\e){-df{e)lde)de 

where A is a constant, t c the hopping probability in the c-axis direction, N(e) 
the DOS, and f(e) the Fermi function. For the spin susceptibility, we assume 
simply the Pauli paramagnetism 

Xpauu = A P_Jf{e)f^-)de (8) 

where \xb is the Bohr magneton. Since the magnitude of Pauli paramag- 
netism also reflects the DOS, the coincidence of T* and T* c simply indicates 
that the pseudogap opens at these temperatures (line "B" in Fig. 1). We 
have confirmed this pseudogap opening at T* /T* c by scanning tunneling spec- 
troscopy (STS) 00] (Some of the data is shown as T^ unnel in Fig. 18). This 
result contradicts the explanation where the gradual decrease in x is a t" 
tributed to the development of AF-spin correlations, and the long-debated 
origin of the semiconductive p c (T) is now clear: it is the pseudogap opening. 

Our results here may conflict with NMR measurement (Hj and the ARPES 
[3], which show that the pseudogap opens at around T* (line "A" in Fig. 
1). Their characteristic temperatures, (T£ and T\ RPES ) ) are also shown in 
Fig. 18. We could not observe any discontinuous change in temperature 
dependence or anisotropy of the susceptibility at around line "A". We point 
out, however, that if T* in [S] is defined as the temperature at which \jT{T 
departs from high-temperature Curie- Weiss behavior, it will come close to 
our pseudogap phase boundary "B". The leading edge analysis of ARPES 
may be too simple to detect the pseudogap opening temperature precisely, 
since the spectra broadens a lot due to both a high temperature measurement 
and a large pseudogap energy. 

The p c and x ma Y be particularly sensitive to the onset of the pseudogap 
opening. A saddle point is observed for many high-T c superconductors in the 
band dispersion at (vr,0) [JT]. The band near (tt,0) is so flat that it creates 
a large peak [van Hove singularity (vHs)] in the DOS. Therefore, the DOS 
is concentrated near the (tt,0) region. In addition, a recent ARPES study 
has shown that the pseudogap starts to open from the (tt,0) direction [4*2*] . 
These situations explain why p c and x are sensitive to the pseudogap opening. 
Further, a recent band calculation has shown that the hopping probability 
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t c in the c-axis direction depends on momentum as t c oc (cosk x a — cosk y a) 2 : 
t c takes maximum at (vr,0). This momentum-dependent t c will enhance the 
sensitivity of p c to the pseudogap opening. On the other hand, p a may 
be less sensitive to the pseudogap opening. The carriers around (vr,0) are 
exposed to strong scattering ("hot spot"), while those around (7r/2,7r/2) are 
not scattered very much ("cold spot"). Consequently, carriers at the "cold 
spot" mainly contribute to p a . When the pseudogap opens from the "hot 
spot", p a may not practically be affected. p a begins to feel the pseudogap 
opening when it reaches the "cold spot". However, the decrease in the p a 
means that the pseudogap is not the type that simply reduces Fermi surface 
area, such as charge-density-wave formation, but the type that only reduces 
the scattering rate, such as opening a gap in the spin excitation spectrum. 

The experimentally obtained pseudogap phase boundary (T* T* c and 
Ttunnei m Fig. 18, "B" in Fig. 1) is not smoothly connected to the T c 
boundary of the heavily overdoped state; rather, it seems to cross the T c 
boundary at a slightly overdoped level and approach the quantum critical 
point at around p=0.20 (T=0) (HUEI]. This suggests that the pseudogap is 
not simply a precursor of the superconducting gap [3T] , but it is a manifesta- 
tion of some quasi-order that competes with superconductivity [131 EI] . Re- 
cently, several spectroscopic studies heve provided evidence supporting this 
view [321 133 EB E3 ESI H3 HH1 EH] • Tunneling spectroscopy by using grain 
boundary junctions [UJ showed, for the electron-doped Pr 2 _ a; Ce x Cu04_j / and 
La 2 ^ 1 Ce :r Cu04_ ?/ , that a pseudogap state appears for the ground state when 
their superconducting state is destroyed by applying a magnetic field. This is 
important in that they established the pseudogap phase boundary T* within 
the superconducting region. Similar results have been reported for hole- 
doped Bi-2212 by using ITS 031111. In the YBa 2 Cu 3 6+( 5 case, T* c |7J and 
T* [5| indicate that the pseudogap vanishes at the optimum doping level, 
suggesting the boundary is near "A". However, even in this case, the true 
boundary sits at somewhat higher temperatures than line "A" of Fig. 1 and 
eventually crosses the T c boundary at the optimum doping. In this sense, 
YBa 2 Cu 3 6+< 5 is not exceptional. We speculate that the cause for somewhat 
lower pseudogap phase boundary of YBa 2 Cu 3 6+ 5 is the 3D nature specific 
to this system [50j . 

3.1.5 TEMPERATURE/DOPING DEPENDENCE OF X 

In order to analyze the T dependence of anisotropic magnetic susceptibili- 
ties quantitatively, we plotted Xc(T) vs. XabiT) with an implicit parameter 
T for several doping levels (Fig. 19) [27]. In the analysis, we assume that 
the T dependence comes from spin susceptibilities of the doped Cu0 2 plane 
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[eq. (4)]. Each of the Xc(T) vs. XabiT) plots [Figs. 19(a), (b), and (c)] 
shows the same linear relation, i. e., Xc oc l.Gxab, except near T c (<140 K). 
This implies that the spin susceptibility has a single component with the 
ratio of anisotropic ^-factors, (g c /g a ) 2 =l-Q. The Curie-like explanation for 
the x °f overdoped samples may therefore be excluded because it is unlikely 
that the localized paramagnetic centers, which could have appeared addi- 
tively in the overdoped state, would accidentally have the same anisotropy 
as the main spin component. This single-component hypothesis is consistent 
with NMR experiments (SB E2] and ensures that we are actually measur- 
ing the intrinsic spin susceptibilities. The anisotropy ratio coincides with 
that for La 2 _ :E Sr :E Cu04 [53], but is somewhat larger than Johnston's result 
for Bi-2212 $R and the result for YBa 2 Cu 3 6+5 [HH E21 We believe our 
anisotropy ratio in the Bi-2212 is more accurate than in [HU], since our value 
was directly obtained by the Xc{T) vs. XabiT) plot whereas, in [30J, it was 
estimated assuming a simple model. The nonexistence of the CuO chain 
structure may explain the difference between Bi-2212 and YBa 2 Cu306+<5. In 
Figs. 18 and 19, we included the temperatures at which the deviation from 
the linear scaling becomes evident, T sc f. Since the diamagnetic supercon- 
ductive fluctuation effect mainly appears on Xc(T) in the 2D system, T sc f is 
considered to be the onset temperature of superconductive fluctuation. In- 
deed, the diamagnetic component of the superconductive fluctuation Sxdia 
follows the power law, Sxdia oc e~ 2 - 3 , where e = ln(T/T c ), and gives the same 
magnitude independent of the doping level |IU]. This observation justifies 
our interpretation of T sc f. It should be noted that the superconductive fluc- 
tuation is independent of the pseudogap behavior and actually additive in 
the pseudogap region (Figs. 18 and 19). 

As already mentioned, we consider that the T dependence of x f° r under- 
doped samples is caused by the pseudogap opening, which has been evidenced 
by STS |H21 HOI EH IHU - Then, it is natural to consider that the Curie-like 
behavior of x f° r overdoped samples is also caused by a DOS effect, since 
the system is expected to be more Fermi-liquid- like upon doping. One such 
possibility is van Hove singularity (vHs), whose existence has been shown by 
ARPES at somewhat below the Fermi energy (ep) for a slightly overdoped 
sample [41 J. The vHs grows and approaches the ep with increasing doping. 
Such doping-dependent vHs has actually been observed by our STS [471] . 
Figure 20 shows the tunneling DOS taken at slightly above T c for various 
5's. In the heavily overdoped sample with 5 = 0.3, the DOS shows a clear 
peak (vHs) at the Fermi level. Then the peak shifts toward lower energy and 
rapidly broadens with decreasing doping level. For the sample with 5 < 0.27, 
the pseudogap develops as a dip structure at the Fermi level. The dip struc- 
ture becomes more pronounced when the doping level is decreased. Here, 
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we assumed simply the Pauli paramagnetism Xpauii for the single-component 
spin susceptibility [eq. (8)]. The Xpauii reflects the thermally averaged DOS 
near Ep. Therefore, the anomalous energy- dependent DOS (pseudogap and 
vHs) and its doping evolution (Fig. 20) with the assumption of Xpauii will 
explain the temperature/doping dependence of x observed. This kind of 
interpretation was first proposed by Loram et al. |3Tj and later asserted 
by Williams et al. |55j . Our data here support their interpretation from a 
different point of view. 

The assumption of the Pauli paramagnetism for the spin susceptibility 
enables us to estimate the pseudogap value |2Zj- Here in the underdoped 
state, the energy independent background DOS, N , was assumed (The vHs 
sits far below Ep and is smeared out). For the pseudogap, the BCS-type DOS 
(Dynes formula) with the d-wave symmetry gap was assumed: 



where V is the carrier scattering rate, assumed here to be 2fc#T in a consis- 
tent way with the transport data, and A is the maximum d-wave gap, the 
temperature dependence of which was neglected. The isotropic component 
of the spin susceptibilities, x l sp°im was obtained by averaging the anisotropic 
spin susceptibilities, i.e., X i s S p°n=( 2 /3)xtpin+( 1 /3)Xs P in, where xf pin and x%in 
were obtained using the values in the literature for x\b an d Xc EOj- Then 
we can fit the experimental spin susceptibility Xspin(^) D Y the calculated 
ones with a free parameter A as shown in Fig. 21. The pseudogap value A 
was estimated as 65±3 and 48±3 meV for the underdoped (5=0.22) and opti- 
mum doped (5=0.25) samples, respectively. These pseudogap values coincide 
with our direct STS measurements jlOj and, possibly, with the "high energy 
gap" suggested by photoemission spectroscopy (SU EZJ , but are larger than 
the superconducting gap (20 meV<A<45 meV) obtained by ARPES/STS 
|3 EH 113 02] . These different gap values between the pseudogap and the 
superconducting gap suggest again that the two gaps are different in origin. 



We controlled the oxygen content S by annealing a sample while varying Ar 
and O2 gas flow ratio and/or temperatures. The annealing conditions for 



Bi-2223 samples used are (a) 2 5xl0- 3 torr, 600°C; (b) 2 0.01%, 600°C; 
(c) 2 0.1%, 600°C; (d) 2 1%, 600°C; (e) 2 10%, 600°C; (/) 2 600°C; (g) 
2 500°C; (h) 2 400°C; and (i) HIP 2 400 atm, 500°C (These descriptions 




(9) 



3.2 Bi 2 Sr2Ca2Cu 3 Oio +( 5:Bi-2223 



3.2.1 DOPING DEPENDENCE OF T c 



c 
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are used in all figures). The superconducting transition temperature T c was 
defined as the onset of the Meissner effect. 

Figure 22(a) shows the normalized transition temperature T c plotted 
against the relative change of the c-axis length from that of the optimum 
doped sample "/" jHB|- The relation between the T c and c-axis length of Bi- 
2212 is also plotted for comparison. The T c and c-axis length are 89 K and 
30.864 A for optimally doped Bi-2212, and 108 K and 37.119 A for Bi-2223 
(sample /), respectively. The c-axis length monotonically decreases with in- 
creasing 5 both in Bi-2212 and Bi-2223, indicating that oxygen is actually 
incorporated into crystals. In the case of Bi-2212, reflecting the bell-shaped 
doping dependence of T c , which is common behavior in monolayer or bilayer 
cuprates, T c increases with decreasing c-axis length, reaches its maximum, 
and then decreases with decreasing c-axis length. In the case of Bi-2223, 
we see a similar increase of T c with decreasing c-axis length. However, T c 
keeps its maximum value as the c-axis length further decreases. This kind of 
anomalous doping dependent T c was first reported in the four-layered com- 
pound (Cui_ I _5C :c )Ba2Ca3Cu40 y (x=0-0.3) jST?] . 

To confirm the carrier doping in the constant-T c region, we measured the 
doping dependence of the thermopowers S [Fig. 22(b)] (SHj- The magnitude 
of the thermopower monotonically decreases with increasing 5. This result 
clearly shows that the carrier was doped continuously even in the constant-T c 
region. The room-temperature thermopower is considered to be a universal 
measure of the doping level [001 ■ By using this measure, the (average) doping 
level of sample / can be considered the "optimal" doping. Then, the sample 
i would be slightly overdoped with carrier concentration of about p = 0.185. 
Thus, it would show T c of 102 K, if the carriers were homogeneously doped. 
Hereafter, we will refer to the constant-T c region as the overdoped region. 

3.2.2 p ab and p c 

The temperature dependence of in-plane resistivity p a b(T) with various <5's is 
shown in Fig. 23 j^B]- The overall doping dependent behavior, such as the 
absolute values of p a b and the slopes dp a b/dT, is similar to Bi-2212 [26J except 
for T c . T c determined by zero resistivity agrees with those determined by 
magnetic susceptibility (inset of Fig. 23). A difference appears in the residual 
resistivity. At all doping levels, the samples show negative residual resistivity 
as indicated by the solid lines in Fig. 23. High-T c materials with T c larger 
than 100 K tend to show negative residual resistivity, although we do not 
know its relevance to T c . The Bi-2223 also seems to belong to this class. The 
underdoped samples (denoted by b, d, e and /) show a downward deviation 
from high-temperature T-linear behavior below a certain temperature T* b , 
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similarly to that of Bi-2212 [211] • T* b increases with decreasing doping as 
168, 192, 203, and 213 K for the samples /, e, d, and b, respectively, as 
indicated by arrows in Fig. 23. 

Figure 24 shows the c-axis resistivity p c {T) for various £'s [SB]- Similar 
to p a b, we can see that T c is pinned at the maximum value in the overdoped 
region. The overall magnitude of p c decreases with increasing 5, which is 
similar to Bi-2212 fZG\ I27j. Here, we can assume the tunneling model [eq. 
(7)] for p c {T) of Bi-2223 also. Then, the decrease in the absolute value of 
p c with 5 would imply an increase in the in-plane DOS and semiconductive 
behavior would be attributed to the decrease of DOS caused by the pseudogap 
formation. The underdoped samples from a to / show semiconductive p c in 
all temperature regions measured, indicating that the pseudogap opens above 
room temperature. A difference appears in the overdoped region. We can 
see that the T* c remains unchanged (« 220 K) for doping levels higher than 
g, while the absolute value of p c continues to decrease. 

3.2.3 INEQUIVALENT HOLE DOPING AND CHARGE CON- 
FINEMENT 

Figure 25 shows a schematic of the basic structure of the Bi-2223 phase [22] • 
As can be seen in this figure, the Bi-2223 system has crystallographically 
inequivalent Cu02 planes: outer Cu02 planes with pyramidal (five) oxygen 
coordination and an inner Cu02 plane with square (four) oxygen coordi- 
nation. Then, there arises the possibility of inhomogeneous doping among 
inequivalent Cu0 2 planes. On the other hand, in high-T c superconductors, 
there is a consensus that the sets of Cu0 2 planes separated by the blocking 
layer are only weakly coupled and the interaction between them can be un- 
derstood as a tunneling process. This is known as the charge confinement 
effect from the theoretical point of view |61j . However, it is still unknown 
whether the confinement works for the Cu02 planes within the unit cell. 

As shown in the previous section, we have observed in Bi-2223 that T c and 
characteristic temperature T* do not change in the overdoped region, while 
the magnitude of p c monotonically decreases and shows about half the value 
of an optimally doped sample (/) in the oxygen-rich state (i). The anomalous 
out-of-plane transport behavior indicates that there are more holes in the 
outer planes than in the inner plane. That is, the behavior can be understood 
by assuming that the holes are mostly doped in the outer Cu02 planes and 
that the inner Cu02 plane remains at the optimum doping level when the 
doping is increased. The decrease in p c can be attributed to the increased 
electrical coupling between the outer planes. It would be difficult to explain 
such a large decrease in p c if there were more holes in the inner plane than 



19 



in the outer planes, since the p c would probably be dominated by the tunnel 
conductance between the outer planes separated by the blocking Bi 2 2 layer. 
Bulk T c is determined by the always optimally doped inner plane, which has 
higher T c than the outer planes. In this case, bulk superconductivity is 
driven by a proximity effect. This inequivalent hole doping agrees with that 
observed in a four-layered compound (Cuo.6Co.4)Ba2Ca 3 Cu40i 2+ j / (T c = 117 
K) • Our result shows that this inequivalent hole doping nature is generic 
in multilayered high-T c cuprates that have three or more Cu0 2 planes in the 
unit cell. 

On the other hand, T* c pinning in the overdoped region, as well as semi- 
conductive p c , is also determined by the always optimally doped inner plane. 
Since the outer planes are overdoped, they do not solely show the pseudogap 
behavior. Then the observed pseudogap effect may come from the tunneling 
through inner planes. This strongly suggests that carriers in the normal state 
are confined to individual Cu0 2 planes. 

3.2.4 WHY T c IS HIGHER IN TRIPLE-LAYERED COMPOUNDS 

It is well established in mono- or bilayer high-T c superconductors that T c 
first increases, takes a maximum, and then decreases with increasing doping 
p (Fig. 1). Although this dome-shaped relation between T c and p is not 
completely understood, the model proposed by Emery and Kivelson may 
be promising. In this model, the trade-off between two basic ingredients 
of superconductivity, i.e. pairing potential A and phase stiffness p s , are 
considered. In conventional superconductors, A is much smaller than p s , 
and thus A always limits T c . However, in high-T c cuprates, the two energy 
scales are comparable. The A is maximal in the underdoped region and falls 
rapidly with increasing doping, whereas p s increases with increasing doping. 
Optimal T c is obtained at a crossover from a phase-ordering-dominated region 
(underdoped) to a pairing dominated region (overdoped). Extending this 
idea, Kivelson recently proposed a way to raise T c of high-T c superconductors 
|64j . His idea is that higher T c 's can be obtained in an array of coupled planes 
with different doped hole concentrations, such that a high pairing energy is 
derived from the optimally or underdoped plane and a large phase stiffness 
from the overdoped ones. We observed that, in the overdoped region of 
Bi-2223, outer planes are overdoped while an inner plane remains at near 
optimum doped (inequivalent hole doping). Here, the coupled system that 
Kivelsom assumed seems to be realized. We consider that Kivelson's scenario 
could possibly explain the higher T c of Bi-2223 [58J. 

Recently, many spectroscopic measurements have been performed on the 
Bi-2223 single crystals j23 ESJ E33 E3 E3 EH- Fen et al. observed 
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by ARPES that the superconducting gap magnitude A and the relative 
coherence-peak intensity, which is regarded as a measure of the superfluid 
density, both scale linearly with T c for the optimally doped Bi family (Bi- 
2201, Bi-2212 and Bi-2223). Based on these observations, they argued that 
the higher T c of Bi-2223 is caused by the enhancement of both the pair- 
ing strength and phase stiffness, which is consistent with the original idea 
of Emery and Kivelson jHS] • By using short-pulse interlayer tunneling spec- 
troscopy (ITS), Yamada et al. [65J observed, for Bi-2223, a similar or a slight 
larger superconducting gap than that of Bi-2212, although they are skeptical 
about the scaling between A and T c . They also reported doping dependence 
of A, where A monotonically decreases upon doping even in the overdoped 
(T c constant) region, and they did not observe two gap structure in those re- 
gion. This behavior is resonably explained in terms of McMillan's tunneling 
model for the proximity effect provided one assumes a somewhat strong 
coupling between inner and outer planes. On the other hand, Sato et al. [66J 
emphasized the scaling behavior. They showed that the characteristic ener- 
gies not only of superconducting but also of pseudogap behavior are scaled 
with T c among the Bi family. From these observations, they implied that the 
effective superexchange interaction J e ff of Bi-2223 is larger than that of the 
others (Bi-2212 or Bi-2201), consistently with spin mechanism for supercon- 
ductivity. The research on the physics of triple-layered or more generally of 
multi-layered (the number of CuC>2 plane n > 3) high-T c superconductors, 
has just been launched. At present stage, the central issue of why T c is higher 
in those systems is an open question. 

4 SUMMARY 

Using the I-TSFZ method, we have succeeded in growing large high-quality 
Bi-2212 and Bi-2223 single crystals. This method makes the best use of 
the advantage of the TSFZ method that crystal growth can be performed 
at one point on the temperature-composition phase diagram. The I-TSFZ 
method will be especially useful for materials that are difficult to grow. By 
the systematic transport and susceptibility measurements for those Bi-2212 
and Bi-2223 single crystals, we found that (i) the pseudogap phase boundary 
crosses the T c boundary near optimal doping and approaches the quantum 
critical point (Fig. 26), suggesting that the pseudogap is not simply a precur- 
sor pairing, but a manifestation of some (hidden) quasi-order that competes 
or at least coexists with superconductivity; (ii) the semiconductive p c {T) 
is caused by the pseudogap opening; (iii) temperature/doping dependence 
of magnetic uniform susceptibility x can be interpreted by the anomalous 
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energy- dependent DOS [pseudogap and van Hove singularity (vHs)]; and (iv) 
charge confinement occurs at each individual Cu0 2 plane. Furthermore, sys- 
tematic transport measurements for Bi-2223 revealed that T c and T* c become 
pinned when doping is increased in the overdoped region. This anomalous 
behavior is interpreted in terms of the inequivalent hole doping (an always 
optimal inner plane with overdoped outer planes) and the charge confine- 
ment. The combination of a largr superconducting gap A in the inner plane 
and a large superfluid density of the outer planes would explain the higher 
T c of Bi-2223 (Kivelson's model). Our finding suggests that both high T c and 
large phase stiffness, which measures the ability to carry a supercurrent, are 
compatible in multilayered high-T c cuprates. From this point of view, the 
multilayered system would be very promising for practical use. 
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Figure captions 

Fig. I. Schematic phase diagram of high-T c superconductors (for introduc- 
tion). 

Fig. 2. Principle of the flux method, (a) Phase diagram and (b) practical 
image. 

Fig. 3. Principle of the TSFZ method, (a) Phase diagram and (b) practical 
image. 

Fig. 4. Schematic phase diagram of Bi family (Bi-2201, Bi-2212 and Bi- 
2223). 

Fig. 5. (a) Solute concentration n sn ahead of crystal growing in solution, (b) 
Equilibrium liquidus temperature T L ahead of crystal growing in solution. 
Temperature gradient smaller than "B" (for example "A") causes constitu- 
tional supercooling. 

Fig. 6. (a) Obtained Bi-2212 boule. (b) Optical polarization micrograph 
of a transverse section of the grown boule. Large single crystal grain with 
thickness larger than 1 mm is seen in the central region (from bottom left to 
upper right region) of the micrograph, (c) Grown crystal cleaved from the 
boule. 

Fig. 7. Crystal growth of Bi-2212 with a growth rate of 0.05 mm/h. (a) Op- 
tical polarization micrograph of a transverse section of the grown boule. Left 
part is a single domain; right part is heavily multi-domain, (b) Sketch of the 
unstable molten zone, which is likely to drop down, (c) Optical polarization 
micrograph of a solid-liquid interface, (d) Enlarged solid-liquid interface of 
the single domain part in (c). (e) Enlarged solid-liquid interface of the multi 
domain part in (c). 

Fig. 8. Crystal growth of Bi-2223. (a) Premelted feed rod at 25 mm/h. (b) 
Actual crystal growth images near molten zone, taken once every week, (c) 
Obtained Bi-2223 boule. Large crystals were found at the last 1-cm part, (d) 
Optical poralization micrograph of a transverse section of the last 1-cm part 
of the boule. (e) Grown crystals cleaved from the upper right side of (d). 

Fig. 9. DC susceptibility of the as-grown Bi-2223 single crystal. 
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Fig. 10. Equilibrium oxygen partial pressures P02 versus reciprocal temper- 
ature for a Bi 2 Sr 2 CaCu208+,5 single crystal for various oxygen contents (5). 
Triangles, squares, circles, and crosses are experimentally obtained data for 
5 = 0.29, 5 = 0.28, 5 = 0.27, and 5 = 0.26, respectively. 

Fig. 11. DC susceptibilities of the Bi 2 Sr 2 CaCu20 8+( 5 single crystals for 5 = 
0.22, 5 = 0.24, 5 = 0.26, and as-grown. Data were normalized at the lowest 
temperature. 

Fig. 12. Oxygen content dependence of T c in Bi-2212. 
Fig. 13. Furnace for the precise oxygen control. 

Fig. 14. In-plane resistivities p a of Bi2Sr2CaCu20s+5 single crystals versus 
temperature for various oxygen contents (5). The solid straight lines, which 
are linear extraporation of p a at higher temperatures, are eye guides for 
the near optimally doped (5 = 0.24) and underdoped (5 = 0.22, 0.217, and 
0.2135) samples. The temperatures T* a at which the p a deviates from T- 
linear behavior are shown by arrows for the underdoped samples. 

Fig. 15. (a) Out-of-plane resistivities p c (T) of a Bi 2 Sr2CaCu20 8+ 5 single 
crystal (sample I) for various oxygen contents (0.2135 < 5 < 0.28). (b) 
Out-of-plane resistivity p c (T) of a Bi 2 Sr2CaCu208+<5 single crystal (sample 
II) for oxygen content 0.26<<5<0.3. The solid straight lines, which are linear 
extraporations of p c at higher temperatures, are eye guides for the overdoped 
(5=0.26, 0.27, and 0.28) samples. Arrows indicate the temperatures T* c 
below which p c shows a typical upturn. 

Fig. 16. Magnetic susceptibilities (a) Xc(T) (H \\ c) and (b) Xab{T) (H _L 
c) of the Bi 2 Sr2CaCu20s+5 single crystals with various oxygen contents (5) 
under H=5 T. The solid straight lines in (a), which are linear extraporations 
of Xc at higher temperatures, are eye guides for the overdoped (5=0.26, 0.27, 
and 0.28) samples. Arrows indicate the temperatures T* at which Xc starts 
to decrease from its linear high-temperature behavior. 

Fig. 17. Scaling behavior for the resistivity ratio r = p c , P eak/ Pc(RT) versus 
doping levels (p). 

Fig. 18. Characteristic temperatures obtained by various measurements ver- 
sus carrier concentration p. The p values (including for the previous results) 



28 



were systematically estimated from the empirical T c vs. p relation in the 
text. 

Fig. 19. Xc(T) versus Xab{T) plots for (a) underdoped (5=0.22), (b) slightly 
overdoped (5=0.26) and (c) heavily overdoped (5~0.3) Bi 2 Sr 2 CaCu 2 0s+,5 sin- 
gle crystals. Several typical temperatures are shown by arrows. 

Fig. 20. Doping dependence of the tunneling DOS in the normal state. 

Fig. 21. Isotropic component of the spin susceptibilities, Xs»1nCO> °f a 
Bi 2 Sr 2 CaCu 2 8 +<5 single crystal for oxygen content 5=0.22 (underdoped) and 
0.25 (optimum doped). The solid lines are numerical fits, assuming the BCS- 
type DOS for the pseudogap. 

Fig. 22. (a) Normalized transition temperature T c plotted against the c- 
axis variation. The T c and c-axis length of optimum-doping samples are 
89 K, 30.864 A for Bi-2212 and 108 K, 37.119 A for Bi-2223, respectively, 
(b) Thermopower S of Bi 2 Sr 2 Ca 2 Cu 3 Oio+,5 single crystal annealed in various 
atmospheres. 

Fig. 23. In-plane resistivity p a b of Bi 2 Sr 2 Ca 2 Cu30io+«5 single crystal annealed 
in various atmospheres. The solid straight lines, which are linear extrapola- 
tions of p ab at higher temperatures, are eye guides. The temperatures T* ab 
at which the p a b deviates from T-linear behavior are shown by arrows. The 
scale is expanded in the inset for a better view around T c . 

Fig. 24. Out-of-plane resistivity p c of Bi 2 Sr 2 Ca 2 Cu 3 Oio+«5 single crystal 
annealed in various atmospheres. Overdoped behavior of out-of-plane resis- 
tivity is shown in the inset. The solid straight lines in the inset, which are 
linear extrapolations of p c at higher temperatures are eye guides for the over- 
doped sample. Arrows indicate the temperatures T* c below which p c shows 
a characteristic upturn. 

Fig. 25. Crystal structure of Bi 2 Sr 2 Ca 2 Cu 3 Oi +5. Inner square-planar Cu0 2 
plane and outer pyramidal Cu0 2 planes are shown schematically. Bold lines 
denote a unit cell. 

Fig. 26. Schematic phase diagram of high-T c superconductors (for summary) . 



29 



Figure 1, T. Wat ana be et al. 




Doping level (p) 



Figure 2(a), T. Watanabe et al 



) 



Q+Liquid \ 


Liquid 
\ Liquidus 
r*r\ line 

:uration\\ 




T+Liquid , 
supersal 





Q T(Target crystal) 

Composition 



Figure 2(b), T. W a tana be et al. 




Heating 
element 



Figure 3(a), T. Watanabe et al. 



) 



Q +Li ^ uid \ Liquid 




\ Liquidus 








« # V 




T+Liquid j \ 







Q T(Target crystal) s 

Composition 



Figure 3(b), T. Watanabe et al. 



(b) 



Feed rod 



Solvent 
composition 



Seed 
crystal 



Figure 4, T. Watanabe et al. 




Figure 5(a), T. Watanabe et al. 




Figure 5(b), T. Watanabe et al 







B 



(b) 




Distance (x) 



Figure 6(a), T. Watanabe et al 



III I'll I'll IHII1 H;l 




Figure 6(b), T. Watanabe et al. 




Figure 6(c), T. Watanabe et al 




Figure 7, T. Watanabe et al 




Figure 8, T. Watanabe et al 




Figure9, T. Watanabe et al. 



4 - 



~\ 1 — i — | — i — i 1 1 1 1 — i — | 1 1 1 1 1 1 1 — | 1 — i — r 



H//a 30G 



Field. Cooled 



Zero.Field.Cooled 



_i i i i i_ 



j i i i i_ 



20 40 60 80 100 120 

Temperature (K) 



Figure 1 0, T. Watanabe et ah 



IS 




1.2 1.4 1.6 

1000/T (K~ l ) 



1.8 



Figure 1 K T. Watanabe et ah 



0.2 



H//ab 

q \ Field (30gauss) cooled 



-0.2 - 



-0.4 - 



* -0.6 



-0.8 - 



-1 



-1.2 



5=0.22 

5=0.24 

5=0.26 
as-grown 




20 30 40 50 60 70 SO 90 100 



Temperature (K) 



Figure 12, T. Watanabe et aL 




0.22 0.24 0.26 0.28 



0.32 



Excess oxygen (8) 



Figure 13 , T. Watanabe et al. 



Vent. Heating Element 



Slotted Valve 

\ 



Oil 

Diffusion 
Pump 



Quarts Tube 



Value Control 



Vacuum j ^ 
Meter 

Sample 



MFC 



Personal 
Computer 



Temp. 
Control 



Ar 



02 



Figure 14, T. Watanabe et ah 




Temperature (K) 



Figure 15(a), T. Watanabe et al 




Temperature ( K ) 



Figure 15(b), T. Watanabe et aL 





■ \ 

■¥■ — 

i\ \ 




■ i 


AT*p„ ^5=0,26 : 




; AVp ^8=0.27 


*_!_!. 

■ ■ i 
i i , 

■ i i 


; AT * p -6=0.28 : 


i J. ^jtb 1 i i i i 1 i i i i 1 i i i i 1 i i i i 



i J ^JUf i ■ ■ i i i i i ■ i i ■ ■ i i i ■ i i i 

50 100 150 200 250 300 



Temperature (K) 



Figure 16(a), T. Watanabe et ah 




1L1 i! i'-'I i i i i I i i i i I i i i i I i i i i I i i i i I i i i l 



50 100 150 200 250 300 350 400 

Temperature (K) 



Figure 16(b), T. Watanabe et al. 




Temperature (K) 



Figure 17, T. Watanabe et aL 




Figure IS, T. Watanabe et ah 




Figure 19(b), T, Wa tana be et al. 




Figure 19(c), T. Watanabe et ah 




1 1 — i 1 — i 1 — i 1 1 1 1 — | 1 — | 1 — | 1 — | r 




J I I I I I I I I I I I I I I I I I L 



-400 -200 200 400 
Sample V (mV) 



Hgurc 20. T. Watanabe ct al. 



Figure 21, T. Watanabe et aL 




50 100 150 200 250 300 

Temperature (K) 



Figure 22(a), T. Watanabe et al. 



i : 



%• 0.95 I- 
H 



0.9 : 



0.85 : 



i i i i i i i i i i i i i i i i i i i i 

i h g f 0.25 (a , « 

■ • • ^^0.24 ,a) - 

0.26 



0.27 



o.23 rr d 



6=0.28 



0.22 



r c 



Bi-2223 
Bi-2212 



doping 



J L 



J i i i i I i i i i I i i i i L 



J L 



-0.15 -0.1 -0.05 0.05 

c-axis variation (A) 



0.1 



Figure 22(b), T. Watanabe et ah 




Temperature (K) 



Figure 23, T. Watanabe et al. 




Temperature (K) 



Figure 24, T. Watanabe et al. 




Temperature (K) 



Figure 25, T, Watanabe et aL 




Figure 26, T. Watanabe et al. 



\ 

\ 




General pseudogap 






phase boundary 






(mono- or bi layer 






System) 




YBa 2 Qi 3 6+ X \ 






pseudogap phask 


\ Quantum 




boundary' \ 




\ Critical 


Antil 


^rro -~ 


■-^i Point 


-mag 


netic \ 


, V\ (QCP) 


Insu 


E,tor / 


■ i \ / 




/ Superconductor 





Doping level (p) 



